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The available data on J photoproduction are analyzed in terms of pomeron exchange, two 
gluon exchange and photon-gluon fusion models. Allowing the pomeron-quark interaction to be 
flavour dependent and introducing the soft and hard pomerons it is possible to reproduce the data 
at v^>10 GeV and small \t\. The two gluon exchange calculations indicate strong sensitivity to the 
gluon distribution function. The results obtained with the most modern MRST2001 and DL PDF 
reproduce the forward J/^ photoproduction cross section at v^>10 GeV. The calculations with the 
photon-gluon fusion model and with MRST2001 and DL PDF are also in reasonable agreement with 
the data on the total J/*!/ photoproduction cross section. We allocate the J/*? photoproduction at 
low energies and large \t\ to the mechanism different from pomeron or two gluon exchanges. We 
consider that this might be axial vector trajectory exchange that couples to the axial form factor of 
the nucleon. 
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I. INTRODUCTION 



i.e. at Vs<10 GeV. 
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J /"^ production by real and virtual photons provides 
an effective way to test QCD dynamics at short distances 
and to verify QCD inspired models. The pomeron ex- 
change model is one of the no^nperturbative QCD ap- 
proaches traditionally applied 1, |3j S IE to vec- 
tor meson photoproduction at high energies. The cru- 
cial point of the model is the flavour independence of 
the pomeron-quark interaction, which can be studied 
through systematic analysis of the photoproduction of 
light, strange, charm and beauty quarks. 

The perturbativc QCD models, such as two gluon ex- 
change ^ and photon-gluon fusion 10, JJ^, 0, 0, 
El 111 El, provide a direct access to the gluon distri- 
bution function. The J/\l/ electro and photoproduction 
is an excellent tool to test the gluon parton distribution 
function (PDF). 

J/'^ electroproduction has been systematically stud- 
ied within these different approaches. Most modern 
gluon distribution functions, such as MRST2001 Ell and 
CTEQ6 |[18| became available just recently through the 
analysis of J/^ electroproduction at large . Very 
recently^the data on J/"^ electro production were ana- 
lyzed [19| with the photon-gluon fusion model. 

Here we present a systematic analysis of J/^* produc- 
tion by real photons. We collect old data available at low 
ener gies , 2 1. 22 . 2iL .24] as well as data at high ener- 
gies ^ [23 13"^ from HI and ZEUS. We also include 
very recently published results f2^. 

Our aim is to investigate whether J/^ photoproduc- 
tion can be understood in terms of a certain model, which 
is also applicable to the photoproduction of other vector 
mesons. It is also a crucial line of our study whether the 
QCD or pQCD models are able to explain the mechanism 
of J/^E* photoproduction close to the reaction threshold. 



II. THE DIFFERENCE BETWEEN u AND J/* 
PHOTOPRODUCTION. 



The Lo, (j) and J/^' vector mesons, which are the mix- 
tures of uu, dd, ss and cc quarks states, have the same 
quantum numbers. For ideal mixing, 9v is near 35°, 
(/>=ss and J /'^=cc, while the w meson is built up with 
uu and dd. However, the mixing is not ideal and both 
the 4) and J/'^ vector mesons contain a certain fraction 
of light quarks. 

Systematic Regge theory analysis 0, [s^ of w meson 
photoproduction shows that at large invariant collision 
energies, ■y/s>10 GeV, the reaction is dominated by soft 
pomeron exchanges, while at y^>100 GcV additional 
small contribution might come from the hard pomeron. 
At low energies, \/s<10 GeV, the lu photoproduction is 
dominated by the meson exchanges. The Regge model 
with contributions from tt and /2 trajectories reason- 
ably describes the data as well as the standard meson 
exchange model with tt, 77 and a contributions. Particu- 
larly large contribution from tt meson exchange to m pho- 
toproduction at low energies can be understood through 
the large w-^-7r+7 partial decay width, which dominates 
the LO meson radiative decay mode. 

An absolutely different situation holds for the J/^ me- 
son. Following the Regge theory it is believed that at 
high energies the J /'^ meson is produced through soft 
and hard pomeron exchanges, where the last dominates 
at y^>50 GeV. It is hard to classify what mechanism is 
responsible for J/^ photoproduction at low energies, be- 
low Vs_=10 GeV. The J/-^ mass mj~3.097 GeV is below 
the DD threshold =i3.739 GeV and OZI (aHHHllal 
allowed J/'^-^D+D decay cannot not be realized in vac- 
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uum. Instead there are large amount of OZI suppressed 
decay channels involving stable hadrons, hadronic res- 
onances and radiative decays A. large fraction of 
radiative J/^" decay proceeds through the production of 
niulti mesonic resonant and non resonant states. The 
J/4'^7r+7 and J/5'^r7+7 both constitute a negligible 
part of the total radiative decay modes, thus the single 
meson exchanges contribution to J/^ photoproduction 
are expected to be very small. 

However, it might be that the mesonic or multi mesonic 
exchange contribution to 'y+N^J/'^+N is negligible 
and the soft pomeron exchange dominates already start- 
ing from the reaction threshold. 

The difference between ui and J/^* total photopro- 
duction cross section is illustrated by Fig^ where the 
reduced total photoproduction cross section Ojj'Js is 
shown as a function of photon energy, i?^. The upper 
axis of FigH indicates the invariant collision energy -y/s, 
where s—mj^+2m]^E^ with mpf being the nucleon mass. 

The uj photoproduction data are shown by squares and 
were collected in Ref. 0- The circles show the J/^* data 
taken from Ref. where the results available until 

1985 were reviewed. The part of the J/^E* photoproduc- 
tion data j22l | were obtained with the nuclear target and 
therefore should be considered as an average between the 
proton and neutron. The J/^ photoproduction thresh- 
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FIG. 1: Reduced uj (squares) and J/vl/ (circles and triangles) 
total photoproduction cross section as a function of photon 
energy. The uj photoproduction data are collected in Ref. Q. 
The J/\[' data available before 1985 are shown by circles and 
taken from Ref. i22il . The J/^ photoproduction data at high 
energies are shown by triangles and are taken from Refs.|24L 
[2^ . |2b . l2g| . The dashed lines show the calculations with 
soft pomeron exchange alone, while the solid lines are the 
sum of soft and hard pomeron contributions. The upper axis 
indicates the invariant collision energy y^. 
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FIG. 2: Reduced uj (squares) and J/'i/ (circles and trian- 
gles) forward photoproduction cross sections as a function of 
photon energy. The ui photoproduction data are collected 
in Ref.;fl]. The circles show J/^ data from Ref. 22], while 
triangles show J/^' photoproduction data from Refs. 27, 29]. 
The dashed lines show the calculations with soft pomeron ex- 
change alone, while the solid lines are the sum of soft and 
hard pomeron contributions. 



old is EjC^S.lb GeV. The Cornell measurement |20j is 
the only result available now for J/\l/ photoproduction 
close to the reaction threshold. We note that the J/^ 
photoproduction data measured at 9.0<ii^.y<11.8 GeV 
at Cornell were incorrectly reviewed in Ref. !22j. The 
original experimental numbers for the differential J/^* 
photoproduction cross section were published pol | in 
the form dcr/dt=(1.01±0.20)exp[(1.25±0.20)t] nb/GeV^. 
The total photoproduction cross section is therefore 
(T=0.48±0.12 nb, taking into account the minimal four 
momentum transfer squared. 

The triangles in Fig^show the J/^* photoproduction 
data, which appeared after 1985 and are pubhshed in 
Refs. 11 m 111 in m It is clear that the en- 

ergy dependence of 'y+p^uj+p and j+N^J/'if+N is 
different at all photon energies. It is known that this 
difference was discovered after the HI and ZEUS mea- 
surements 25, 26, 27, 28, 29] at ^>50 GeV. 

Fig[21 show the u and J/^* reduced differential pho- 
toproduction cross section, da/dt/^/s, extrapolated to 
zero four momentum transfer squared t—0 as a function 
of photon energy Ej and invariant collision energy ^/s. 
The energy dependence of the forward photoproduction 
cross sections again indicates a strong difference between 
the UJ and J/"^ data. Our meaning of the forward photo- 
production cross section is an extrapolation of da/dt to 
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FIG. 3: The slope, b, of the t dependence of the differential 
u (squares) and J/'i' (circles) photoproduction cross sections 
as a function of photon energy, E^, and invariant collision 
energy, The data are from Refs.P, 22, 27, 29]. The 

solid lines show the calculations with soft pomeron exchange 
alone for cu photoproduction and with the sum of soft and 
hard pomeron contributions for J/^ photoproduction. The 
calculations show the energy dependence of the local slope at 
the indicated four momentum transfer squared, t. The dashed 
lines show the fits using a constant value. The J/'i/ photopro- 
duction data were fitted separately at photon energies below 
and above 100 GeV. 



t=0. 

Another value of interest is the exponential slope b of 
t dependence of the differential photoproduction cross 
section, which is defined as 



J/^ meson photoproduction might indicate a rather dif- 
ferent mechanism for their production. Although, the 
quantum numbers of uj and J/^P are the same, they 
have different quark structures. The difference between 
the Lo and J/^t photoproduction data can be attributed 
to the difference between the properties of light and 
charm quarks, such as their interaction through the 
gluon exchange, confinement and fusion into real hadrons 
and photon quark-antiquark fluctuation. Therefore, cur- 
rently available theoretical models for u; and J/^P meson 
photoproduction are inspired by QCD. 

In the next section we review some of the models that 
recently were successfully applied to the J/^* electropro- 
duction off the nucleon. We will compare these models 
with data on J/^E* photoproduction by real photons. Fur- 
thermore, we consider all available data on J/^> photo- 
production in order to make our study in a more system- 
atic way. 



III. POMERON EXCHANGE MODEL. 

The basis of the model is the factorization of the exclu- 
sive vector meson photoproduction amplitude in terms of 
the product of 'y^q+q fluctuation, the scattering of the 
qq system by the proton and finally the qq hadroniza- 
tion into a vector meson. The interaction between the 
qq state confined in the vector meson and the nucleon is 
taken to be given by pomeron exchange. The pomeron 
propagator G{s,t) is parameterized as 



G{s,t)=^[^y^'^ exp(-f KO-1]), 



(2) 



and therefore the pomeron trajectory a(t) drives the ^/s 
energy dependence of the reaction amplitude. The evolu- 
tion of t dependence with energy is controlled by pomeron 
exchange trajectory as well. 

The interaction between the pomeron and quark is con- 
sidered to be a flavour independent. The flavour indepen- 
dence of the quark-pomeron interaction might provide a 
unique way to construct a unified model for vector meson 
photoproduction. 



A. Soft pomeron. 
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X exp(6i). 



(1) 



and is shown in Fig|31 The lu photoproduction data can 
be well fitted with a constant slope fe=7.6±0.1 GeV"'^ 
over a wide range of available photon energies. In 
contrast, the J/^ photoproduction data show a strong 
energy dependence of the exponential slope. At 
Ej<100 GeV the data can be well reproduced with a 
constant slope 6=2.6±0.4 GeV~^, while at high energies 
&=5.1±0.2 GeV~2 

Finally, the available experimental results on uj and 



The pomeron exchange amplitude T for vector meson 
photoproduction is depicted in Fig^^) and is explicitly 
given H H H H H in the form 



r ^3iFi(t) —- — —^{e-ev) 



4m2 



xexp(-y[a(i)-l]) 



Q(t)-1 



4m^ 



(3) 



where Cq and are the charge and mass of the quark, 
e and ey are the polarization vectors of the photon and 
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vector meson, respectively and fv is the meson decay 



constant given by V—i-e 
as 



radiative decay width Fg 



3mi 



(4) 



where a is the electromagnetic coupling constant and my 
is the vector meson mass. 

Furthermore, t is the squared four momentum transfer 
and Fi (t) is the proton isoscalar EM form factor given as 



Flit) 



Ami - 2.8t 



1 



4m2 - t 



(5) 



with irip being the proton mass and to=0.71 GeV^. The 
Fi(t) was introduced assuming the pomeron resembles 
an isoscalar photon. 

The pomeron trajectory a{t) for soft pomeron ex- 
change is given as 



a{t) = 1.08 + a'i, 



(6) 



with a' =0.25 GeV^^. The constant so is not well defined 
theoretically and can be taken as SQ=l/a' following the 
dual model prescription [s^ . 

The free model parameter, determines the strength 
of the effective pomeron coupling to the quark, while fiq 
accounts for the possibility that the coupling to an off- 
shell quark is not pointlike but dressed with the form 
factor given by the last term of Eq.©. Both Pq and, in 
principle, fiq might depend on the flavour of the quark 
and can be fixed by experimental data. 

The phase of the soft pomeron exchange amplitude 
was not written explicitly in Refs. Q, |^ m> Q II] since 
no interference with other exchanges was considered. 
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FIG. 4: The diagrams for pomeron exchange (a), two gluon 
exchange c,d) and photon-gluon fusion model b). 



The amplitude is normalized so that da/dt=a\T\'^. Fi- 
nally, the differential 'f+p^V+p cross section due to soft 
pomeron exchange is given for a real photon as. 



da 



81 mf. 



fit y^t 



7ra 



So 



2qpj (t)-2 



(7) 



where it was assumed that mv—^mq. 

The coupling constant (3q and the form factor fiq were 
fitted |23] through the systematical analysis of uj meson 
photoproduction as 



Pq = 2.35 GeV" 



1.1 GeV. 



(8) 



It was found that the parameter l3q evaluated from 
uj photoproduction was different from Pq=2.0 GeV~^ 
needed in the description [23, Ull of p photoproduction 
data. This discrepancy could not be explained in a sat- 
isfactory way, since once the parameters of the pomeron 
exchange are fixed by p meson photoproduction, they 
should be applicable to the other vector mesons as well. 
Moreover, at this stage we still do not discuss the flavour 
dependence of the quark-pomeron interaction, because 
we are dealing with p and uj photoproduction, i.e. light 
qq states. 

An ideal illustration of this finding is given by the ratio 
of the ^+p—^uj+p and "f+p—^p+p cross sections, shown 
in Fig. as a function of photon energy E-y and invariant 
collision energy y/s. Following Eq. the ratio of oj and p 
meson photoproduction is given by the ratio of uj—>e~^e~ 
and p^e+e^ decay widths 



R 



0.60±0.02 keV 
6.77±0.32 keV 



0.088±0.005, 



(9) 

which is shown by the solid line in Fig. El and considerably 
underestimates the experimental results. 

The dashed line in Fig. |S1 shows the fit to the data by 
a constant value i?::=0.115±0.003, which is close to the 
SU{3) estimates for the uj and p coupling to the photon. 
Finally the difference between the Regge model and the 
data accounts for a factor of ~ 1.3. It is already clear 
that both p and uj photoproduction data could not be 
reproduced by the soft pomeron exchange simultaneously 
with the same set of parameters (3q and pq. 

Now, the dashed line in Fig. ^ shows the contribution 
from soft pomeron exchange to the total reduced uj me- 
son photoproduction cross section. The calculations were 
done with the parameters given by Eq.(|SJ|. The compar- 
ison to differential j+p^uj+p cross section at different 
photon energies is given in Ref. j30|. The contribution 
from soft pomeron exchange to forward uj meson photo- 
production is shown in Fig. |21by the dashed line. 

Obviously, the soft pomeron exchange alone will repro- 
duce the UJ photoproduction data at \/s>6 GeV. At lower 
energies the dominant contribution to the j+p^uj+p re- 
action comes from meson exchanges. 
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FIG. 5: The ratio of 'y+p^uo+p to cross sections 

as a function of plioton energy, E~f, and invariant collision 
energy, y^. The solid line show the ratio of the a;— >e+e~ and 
p^e+e~ decay widths. The dashed line shows the fit to data 
with a constant value of _R=0.115±0.003. 



The t dependence of the photoproduction cross sec- 
tion arising from the soft pomeron exchange is given by 
the proton EM form factor squared F^{t), the quark- 
pomeron form factor squared and by the pomeron ex- 
change trajectory (s/so)^"''*-'^^. The exponential slope, 
&, of the t dependence can be evaluated as 



dt 



da 
It 



(10) 



from the differential photoproduction cross section given 
by Eq.Q. 

The slope due to the proton isoscalar EM form factor 

is 



bi 



5.6 



4m2-2.8t Aml-t to - 



(11) 



and since it depends on t it can be considered as a local 
slope. At t=0 the slope from proton EM form factor is 
~4.6. The slope due to the quark-pomeron form factor 
squared is given as 



(12) 



t 



-t' 



and depends on the cutoff parameter /x, and mass of 
the produced vector meson my. For oj meson photo- 
production at t=0 the slope associated with the quark- 
pomeron form factor is ~4.0. The energy dependence of 



the slope comes from pomeron trajectory and is given for 
soft pomeron exchange as 



63 = 0.51n(0.25s). 



(13) 



Finally, the total exponential slope, b, of the t depen- 
dence for Lu meson photoproduction at t=0 has a mini- 
mum value of 6=8.6 GeV^^ and increases logarithmically 
with energy. This result is in excellent agreement with 
the data shown in Fig|21 Here we also show the local 
slope b at different squared four momentum transfers t. 
It is also clear that the appearance of the proton EM 
form factor in the reaction amplitude dictates quite a 
large value of the minimal available b at t—0. 

The contribution from soft pomeron exchange to J/^ 
photoproduction is given by Eq. (0) . Taking the parame- 
ters (3q and fiq as determined from the data on to meson 
photoproduction one can easy estimate the ratio of for- 
ward J/"^ and uj photoproduction cross section as 



R 



mj(m?j 



=0.13, (14) 



with rj^e+g-=5.26 keV. The experimental results 
shown in Fig [5] indicate that this ratio is of order 0.005 
at y^ic^lO GeV. Therefore the soft pomeron exchange 
with fixed parameters for the quark-pomeron interaction 
substantially overestimates the data on J/'^ photopro- 
duction. However, since it was already realized that a 
simultaneous description of both p and uj photoproduc- 
tion data requires a different coupling constant Pq, one 
might proceed in a similar way and readjust Pq to the 
J/'^ photoproduction data. 

Considering the flavour dependence of the interaction 
between the pomeron and quark one might distinguish 
the light quark coupling /3„=/3d=2.35 GeV~^, which is 
taken to be the same as evaluated from ld meson photo- 
production, and the pomeron coupling to charm quark 
(3c- In that case the 5^ term in the amplitude given by 
Eq.© should be replaced by product bubc- In addition 
the differential cross section of Eq.Q must also be cor- 
rected. 

We found that the choice of /3c=0.45 GeV~^ provides 
a reasonable description of the data around y/s~\Q GeV. 
The dashed lines in Figs lll2l show our calculation of the 
soft pomeron contribution to the reduced total and for- 
ward J/^ photoproduction cross section. 

However, resolving the problem with the reduction of 
the soft pomeron contribution to J/"^ photoproduction 
leaves one far short of a description of the J/^ meson 
data collected in Figs lll2i;II 

First, although the soft pomeron seems to reproduce 
the total J/^E* photoproduction cross section starting just 
from the reaction threshold, the model absolutely fails to 
describe the forward cross section and slope b of the t 
dependence at low energies. Recall, that the minimal 
value of the slope is given by the proton scalar EM form 
factor and this is already larger than the experimental 
results at low energies. 
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Second, the energy dependence of the total and for- 
ward J/'^ photoproductfon cross section and the slope h 
at high energies cannot be reproduced by soft pomeron 
exchange. 



a'~0 GeV-2. 

The hard pomeron trajectory was introduced in Ref. 



as 



B. Hard pomeron. 

There are different attempts to reproduce the energy 
dependence of the total J/^ photoproduction cross sec- 
tion without introducing a new pomeron trajectory. The 
strategy is to readjust the t dependence in such a way 
that it becomes flatter, providing the correct total cross 
section after integration over the large t region accessi- 
ble at high energies. One of the more natural ways is to 
introduce a trajectory that is nonlinear in t, instead of 
that given by Eq.|(Bl). Obviously the J/^* photoproduc- 
tion data at t=0, shown in Fig[21 do not support such 
solution. 

Instead, the forward J/^* photoproduction cross sec- 
tion requires that at t=Q the cross section is proportional 
to ^/s, because the reduced cross section in Fig|51is al- 
most independent of energy for \/s>10 GeV. From Eq.lj^J 
one can show that the trajectory necessary to reproduce 
the data has an intercept Q;(i=0)~1.25. Moreover, since 
the energy dependence of the slope, 6, is proportional to 
a'ln(s), the data shown in Fig|31 support the value of 

7 + p ^ J/4^ + p 
.3, 



a{t) = 1.44 + O.li. 



(15) 
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FIG. 6: The J/$-|-p differential cross section as a func- 

tion of four momentum transfer squared, t, measured at differ- 
ent invariant collision energies \fs by HI Collaboration [27| . 
The solid lines show the calculations including both soft and 
hard pomeron exchanges. 



The calculations with both hard and soft pomeron ex- 
changes, including their interference, are shown by the 
solid lines in Figs lll2l5l The coupling constant be- 
tween the hard pomeron and charm quark was fitted as 
/3c=0.05 GeV~^. The calculations reproduce quite well 
the data on total and forward J/^* photoproduction cross 
section and slope b oit dependence at ^/s>10 GeV. Note 
that for both soft and hard pomeron exchanges the par- 
tial slope &i, coming from the form factor in the upper 
pomeron-quark vertices, is 6i~0.37, because of the large 
J/'f mass appearing in Ea. Hll|) . It is clear that the hard 
pomeron dominates at high energies. 

FigsElshows the differential cross section for J /'^ pho- 
toproduction at high energies. The data available at 
small t can be described well by the model. 

While the introduction of a hard pomeron exchange 
allows one to resolve inconsistencies at high energies, 
the discrepancies between the model and the data at 
a/s<10 GeV still remain. 

The most crucial data are shown in Figd where we col- 
lect available data on differential ^+N~*J /'^-\-N cross 
section at relatively small energies, -ys<23 GeV. The 
solid lines show the calculations with soft and hard 
pomeron contributions. The dashed lines show the fit 
to the soft part of the differential spectra. 

The comparison between the differential spectra and 
calculations at -y/s<10 GeV illustrates that both slope 
and extrapolated cross section at i=0 could not be re- 
produced by the model. The rough agreement between 
the integrated theoretical and experimental cross sections 
shown in Fig^ should be considered as an accident. The 
pomeron exchange model is far from describing the data 
at Vs<10 GeV. 

The comparison at y^~20 GeV shows that pomeron 
exchange allows one to describe the steep dependence 
of the differential ^~\-N^J J'^+N cross section at small 
t. At large four momentum transfer squared there are 
additional contributions resulting in a soft component of 
the spectrum. It is worthwhile to note the slope of the 
soft component observed at -^8^20 GeV is close to that 
measured at lower energies and is around 0.7-^2.8 GeV~^. 

FiglHlshows recent data on differential j+N^J/^^+N 
cross section collected by ZEUS Collaboration |23| at 
30<^/s<170 GeV. The solid lines show the calculations 
with soft and hard pomeron contributions, which well 
describe the data at low \t\. At large four momentum 
transfer squared we detect a soft component of the spec- 
trum UlEili^. 
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FIG. 7; The 7+p— >J/^'+p differential cross section as a func- 
tion of four momentum transfer squared t measured at differ- 
ent invariant collision energies ^/s. The data are from Cor- 
nell Hal, E687 2^, EMC 23] and SLAC Collaboration l2l . 
The solid lines show the calculations including both soft and 
hard pomeron exchanges. The dashed lines indicate the fit to 
the soft part of the spectra. 



C. Short comment. 

It is clear that assuming the appropriate flavour de- 
pendence of the interaction between pomeron and quark 
the Regge model allows an excellent description of the 
high energy data on vector meson photoproduction. 

There are attempts to recover flavour independence 
by accounting for the mass-dependent corrections, such 
as Fermi motion and quark off-shellness in the qq loop. 
Indeed the vector meson coupling to qq introduced 
through Eq.Q) was obtained within an on-shell approx- 
imation, which is not correct in the case of diffractive 
photoproduction. However, it was found that calcula- 
tions (42i. A3l t44 | with realistic vector meson wave func- 
tions do not allow one to describe simultaneously the p 
and J/^ meson electroproduction, similar to what we 
have found in the present study of uj and J/^* photo- 
production by real photons. It was proposed to use a 
Gaussian wave function instead and to regulate its mean 
value individually for each vector meson. It is not clear 
whether such degree of freedom is better than the indi- 
vidual choice of coupling constant /3g. 



IV. TWO GLUON EXCHANGE MODEL. 

The exclusive vector meson photoproduction ampli- 
tude is again factorized in terms of the 7— fluctu- 
ation, the scattering of the qq system by the proton and 
finally the qq hadronization into the vector meson. The 
qq system interacts with a nucleon through two gluon 
exchange. 

In lowest order perturbative QCD the photoproduction 
amplitude is given by the sum of two diagrams depicted 
in Fig^,c) and can be written as |7i,iS,|9( 

^ = ^ mq as Cqtv F2g (t) 

X jdfDl{l)[D+{l)-D^{l)]G(l), (16) 

where the integration is performed over the gluon trans- 
verse momentum, and are the charge and mass of 
the quark respectively, while as is the QCD coupling 
constant p35| . The meson decay constant is given by 
Eq.Q) The gluon propagator Dg{l) is taken as l/t^, the 
D^{1) is the propagator of the off-shell quark in the di- 
agram where each gluon couples to a different quark of 




0.5 1 1.5 0.5 1 1.5 

-t (GeV^) 

FIG. 8: The 7-l-p^ J/'I'-l-p differential cross section as a func- 
tion of four momentum transfer squared t measured at dif- 
ferent invariant collision energies \/s by ZEUS Collabora- 
tion 29] . The solid lines show the calculations including both 
soft and hard pomeron exchanges. 
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the vector meson and is given as 

D_{1) = (-2to2_2/2)-i. 



(17) 



When the two gluons couple to the same quark the off- 
shell quark propagator is 



-(0 = (-2™^)" 



(18) 



Furthermore, F2g{t) accounts for the t dependence of 
the amplitude given by a two gluon correlation in the 
proton. The form factor F2g{t) is not defined by the 
model and, as proposed in Ref. |3 can be taken to be the 
proton isoscalar EM form factor, -Fi(i), given by Eq.®. 

Function G{1) defines the probability of catching the 
two gluons with momenta I from the proton and it is 
related to the conventional gluon distribution function 
g{x) as nil 



(19) 



Finally the two gluon exchange amplitude becomes 



r 



iV2 7r2 



xgi^: Ql) 



dP 



IqCKs eqfv F2g{t) 

dxg{x, P 



nP^rrP + P) dP 



(20) 



The amplitude is normalized so that da / dt=a\T\'^ and 
in lowest order the J/^t* photoproduction cross section is 
given as 



da 



[xg{a 



(21) 



where x=rrPj/s and Qo~mv- The cross section depends 
on the gluon distribution function squared and the en- 
ergy dependence of da/dt is directly given by the x de- 
pendence of xg{x). 

To compare the two gluon exchange model with data 
on J/^! photoproduction we need to specify the gluon 
parton distribution function. The PDF used in the fol- 
lowing calculations are shown in Fig|3 as a function of 
x and for Q'^—mPj. The upper axis of FigO indicates 
the relevant invariant collision energy, ^—mj/^/x. The 
old GRV98 m and MRSG ^ functions show a very 
steep dependence at small x, while the most modern 
MRST2001 [13 and CTEQ6 ^1^] functions are almost 
flat in X. The DL structure function shows some av- 
erage dependence between the two extreme GRV98 and 
MRST2001 cases. 

The PDF shown in Figl^are very similar for a;>0.01, 
which corresponds to the invariant collision energy 
of^/s<30 GeV. That means that the data available be- 
fore HI and ZEUS measurement are not sensitive to the 
choice of PDF. One might expect that the differences 
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FIG. 9: The gluon distribution functions as a function of x 
at Q^=m?j. Upper axis shows the relevant invariant collision 
ener gy g iven as \/s=mj / \fx. The functions are taken from 
RefslH III in S, . 



between the different gluon functions in reproducing the 
7-l-p— > J/^'-l-p reaction should appear only at high ener- 
gies. Since xg[x^ Qq) enters the photoproduction cross 
section as an unintegrated function squared, the sensi- 
tivity of the energy dependence of J/^t photoproduction 
to the type of PDF should be extremely high. 

The calculations were done for two extreme cases given 
by GRV98 and MRS2001 and for the DL function. The 
energy dependence of the QCD coupling constant ois was 
taken from Ref. 35] . The mass of the charm quark allows 
for additional freedom in the absolute normalization of 
the calculations and we used mg=1.3 GeV. FigllOlshows 
the forward J/5'-|-p cross section as a function of 

photon energy in comparison to the two gluon exchange 
model calculations by Eq. H21|) with different gluon distri- 
bution functions. 

As expected, the different PDF equally well reproduce 
the data around ^s-^Wi GeV. The GRV98 function is 
too steep at small x to describe the J/vP photoproduc- 
tion at high energies. The DL function reasonably well 
reproduces the data over a wide range of photon ener- 
gies, but has incorrect threshold dependence, because of 
its asymptotic behaviour at large a;~l, as can be seen 
from Figl^l 

The modern MRST2001 gluon distribution function 
describes the energy dependence of the forward J/vP pho- 
toproduction cross section very well, providing as well the 
correct threshold behaviour. R is also important to note 
that the two gluon model calculations with MRST2001 
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leave some room for contribution from other processes at 
^/s<8 GeV. As is illustrated by FigsHElandinithe data 
actually require a different photoproduction mechanism 
close to threshold, which certainly can be added to the 
two gluon exchange contribution. 

The t dependence of the two gluon exchange is intro- 
duced artificially through the proton EM form factor Fi 
providing the slope &=4.6 GeV~^ at t=0. Although this 
result is in a agreement with the data shown in FigOone 
may consider this agreement as a reasonable guess. 



V. PHOTON GLUON FUSION MODEL. 

In the photon-gluon fusion model [13 [H El El El 
Esl ll^ (PGF) the photon fuses with a gluon from the 
nucleon and forms a cc pair, as is shown in diagram in 
Fig^). Since a gluon transforms as a color octet, the 
produced cc pair must radiate a soft gluon before the 
hadronization into the color singlet final state. Therefore 
the PGF presumably consider an inclusive charmonium 
or open charm photoproduction associated with excita- 
tion of the final nucleon. 

The model convoluted the the gluon momentum distri- 
bution and the photon-gluon cross section a-yg^cs, which 



IS given as 
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FIG. 10: The forward J/4' photoproduction cross section |22 . 
|2^ . |29| as a function of photon energy, E-^ . The lines show 
the pQCD calculations by two gluon exchange model with 
the gluon distribution functions GRV98 MRST2001 [T^] 
and DL The upper axis indicates the invariant collision 
energy. 



27ra eia 



X In 



2ml)] 

-/3[S2 +4Sm2]), 



1-/3 



(22) 



where Cc and rric are the charge and mass of the charm 
quark, respectively, as is the QCD coupling constant, a is 
electromagnetic coupling constant, s is squared invariant 
mass of the photon-gluon system and 



1 



4m^ 



(23) 



The PGF does not consider explicitly the formation of 
the final state, which should enter through the J/^* wave 
function as in the case of two gluon exchange model. In- 
stead the convolution of photon-gluon cross section with 
the gluon momentum distribution is performed over the 
range from charmonium production threshold, 2mc, up 
to the open charm production threshold 2m d, given by 
the mass of D meson. Finally the charmonium photo- 
production cross section is 



" d~ 

4m^ 



9{x=-,Qo), 



(24) 



where the factor, /, is an adjustable parameter that ac- 
counts for the fraction of the specific charmonium bound 
states available in the mass region between 2mc and 
2to£i. For the 7/4" meson we fit it as /=0.062 with 



=1.3 GeV. Furthermore, we take Qq= 



-m 



The PGF calculations with different gluon distribu- 
tion functions are shown in Fig llll together with data 
on the total j+p^J/^+p reaction cross section. Ob- 
viously, the PGF is less sensitive to the type of gluon 
distribution functions in comparison to the two gluon 
exchange model. The calculations with MRST2001 [l^ 
and DL Ell reproduce the data better than the re- 
sults obtained with GRV98 , which indicates a much 
steeper energy dependence. However the MRST2001 
and GRV98 distributions both have correct threshold be- 
haviour, while DL - not. In principle, the DL gluon dis- 
tribution function can be corrected in order to get the 
correct asymptotic behaviour as x—*l, by multiplying it 
with (1 — x)^. However, then the agreement with the 
data disappears. 

The PGF model does not contain a t dependence and in 
the best case it can be artifically introduced through the 
proton EM form factor. Although the PDF calculations 
are able to reproduce the cross section for J/'i> photopro- 
duction starting from threshold we would not claim that 
pQCD is able to describe the data at •\/s<10 GeV. We 
consider this good agreement between the data and PGF 
calculations due to the correct asymptotic behaviour 
of the gluon distribution function as x— >1, given by 
~{l—x)^ dependence. Figllll shows this function alone 
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in order to illustrate its threshold dependence. Further- 
more, in the two gluon exchange model the PDF enters 
squared and the near threshold behaviour of the J/'^ 
photoproduction cross section is then given by function 
~(1— a;)^*', which is shown in Fig 1 101 However, now the 
dependence is too steep to reproduce the data near the 
threshold. 



VI. CONCLUSION. 

We have analyzed the available data on J/'^ photo- 
production in terms of pomeron exchange, two gluon ex- 
change and photon-gluon fusion models. 

Allowing the pomeron-quark interaction to be flavour 
dependent and introducing the soft and hard pomerons 
it is possible to reproduce almost all available data on 
J/5' photoproduction at ^>10 GeV. At the same time, 
the model reproduce well the data on lo meson photopro- 
duction, again taking into account that the interaction of 
the pomeron with light and charm quarks is very differ- 
ent. However the pomeron exchange model cannot not 
reproduce the data at y/s<lQ GeV and differential spec- 
tra at large The agreement between the calculations 
and total J /"^ photoproduction cross section near thresh- 
old should be considered as accidental, because both the 
slope b of the t dependence and the forward extrapolated 
cross section are in strong systematic disagreement with 
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FIG. 11: The total J/^ photoproduction cross section 
m, m m m as a function of photon energy, E^,. 
The lines show the pQCD calculations by photon-gluon fu- 
sion model with the gluon distribution functions GRV98 |45|l . 
MRST2001 [l3] and DL [l^. The upper axis indicates the 
invariant collision energy. 



the pomeron exchange model at ^<\Q GeV. 

The two gluon exchange pQCD calculations indicate 
strong sensitivity to the gluon distribution function. The 
results obtained with the most modern MRST2001 [l7| 
and DL PDF reproduce the forward J/^* photopro- 
duction cross section at \/s>10 GeV quite well. The 
calculations with MRST2001 show correct threshold de- 
pendence of the cross section and leave room for a contri- 
bution from other mechanisms at low energies. The two 
gluon exchange model does not contain any t dependence. 
Its artificial introduction through the isoscalar EM pro- 
ton form factor allows one to reproduce the slope of the 
t dependence at high energies, but fails at y^<10 GeV. 

The calculations with the photon-gluon fusion model 
again show sensitivity to the gluon distribution function. 
The PGF results calculated with MRST2001 17] and 
DL 19] PDF are in reasonable agreement with the data 
on the total J/^" photoproduction cross section. The 
PGF does not constrain the t dependence. 

Ahhough the PGF calculations with MRST2001 repro- 
duce the threshold behaviour of the total J/^'-fp 
cross section it does not prove that pQCD is able to de- 
scribe the reaction at -\/s<10 GeV. It is just an artifact 
of the correct PDF asymptotic behaviour at x— >1, given 
by ~(1— a;)^ dependence. For instance it is not a case 
for two gluon exchange model, where the PDF enters 
as squared and finally the ~(1— 2:)^" dependence is too 
steep to reproduce the J/^f photoproduction data near 
the threshold. 

We allocate the J/^* photoproduction at low energies 
and large \t\ to the mechanism different from pomeron 
or two gluon exchanges. The energy dependence of the 
differential cross section, da/dt, indicates that this mech- 
anism does not depend on s. Furthermore, small slope b 
of the t dependence indicates that the exchange does not 
couple to the isoscalar EM proton form factor. We con- 
sider that this might be axial vector trajectory exchange 
that couples to the axial form factor of the nucleon. How- 
ever, a definite conclusion requires more detailed study. 
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